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This  work  tries  to  determine  whether  it  is  possible  to  stabilize  C02  emissions  under  a  rapid  increase  of 
Gross  Domestic  Product  in  a  medium  term.  The  paper  is  an  effort  to  study  in  detail  how  changes  in  the 
driving  forces  of  the  economy  affect  C02  emissions.  We  study  the  case  of  Venezuela  for  the  period  1980- 
2025,  using  the  methodology  proposed  by  Robalino-Lopez  et  al.  (2014,  2013)  [1,2],  which  is  based  on  an 
extension  of  the  Kaya  identity  and  on  a  GDP  formation  approach  that  includes  the  effect  of  renewable 
energies.  We  test  the  environmental  Kuznets  curve  (EKC)  [3]  hypothesis  in  a  coming  future  under 
different  economic  scenarios  considering,  therefore,  not  only  past  data  but  also  projections  for  the 
coming  years.  We  use  cointegration  techniques  [4]  and  the  Jaunky  (2011)  [5]  specification  to  test  the 
fulfillment  of  the  EKC  hypothesis  in  Venezuela.  Our  predictions  show  that  Venezuela  does  not  fulfill  the 
EKC  hypothesis,  but  however,  the  country  could  be  on  the  way  to  achieve  environmental  stabilization  in 
the  medium  term.  This  stabilization  could  be  accomplished  combining  economic  growth  with  increasing 
use  of  renewable  energy,  appropriated  changes  in  the  energy  matrix,  and  in  the  productive  sectoral 
structure. 
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1.  Introduction 

There  are  multiple  factors  that  influence  the  level  of  C02 
emissions,  among  them,  economic  development,  population 
growth,  technological  change,  resource  endowments,  institutional 
structures,  transport  models,  lifestyles  or  international  trade  [6,7], 
Latest  observations  show  that  global  C02  emissions  are  far  from 
stabilizing  and  have  experienced  significant  growth  during  the  last 
century.  Clearly,  one  factor  behind  this  increase  is  the  growth  of 
economic  activity.  The  identification  of  the  kind  of  sources  of 
Green  House  Gases  (GHG)  emissions  and  of  its  magnitude  is 
essential  information  for  economic  planning,  for  decision  makers 
and  for  environmental  development.  C02  is  by  far  the  main 
contributor  to  antropogenic  GHG  emissions;  according  to  IPCC 
[8],  C02  represents  76.7%  of  the  GHG  emissions  (approximately 
56.6%  is  from  fossil  fuels,  17.3%  from  deforestation,  and  2.8%  from 
other  sources). 

To  address  this  issue  the  Intergovernmental  Panel  on  Climate 
Change  (IPCC)  has  developed  several  methods  to  estimate  GHG 
emissions,  such  as  the  Reference  Method  [8,  vol.  2,  Chapter  6],  This 
is  a  top-down  technique  that  uses  data  from  the  countiy's  energy 
supply  to  calculate  C02  emissions,  mainly  from  the  burning  of 
fossil  fuels.  This  is  a  straightforward  method  that  can  be  applied 
on  the  basis  of  the  available  energy  supply  statistics  in  most  of  the 
countries. 

The  adoption  of  environmentally  sustainable  technologies, 
improving  energy  efficiency,  saving  energy,  forest  conservation, 
reforestation,  or  water  conservation  are  some  of  the  most  effective 
ways  to  address  the  climate  change  issue  [9],  Furthermore,  Chien 
and  Hu  [10],  Robalino-Lopez  et  al.  [1,2]  show  that  increasing  the 
use  of  renewable  energy  improves  the  C02  emission  efficiency  of 
the  economy. 

The  economic  growth  has  a  direct  impact  on  C02  emissions  of  a 
given  country.  Most  of  the  literature  suggests  that  the  relationship 
between  economic  growth  and  C02  emissions  is  established 
through  the  energy  use  and  influenced  by  the  energy  efficiency, 
the  promotion  of  renewable  energies,  and  by  the  productive 
sectoral  structure.  The  detailed  study  of  the  later  relationship  will 
be  one  of  the  goals  of  this  work.  In  particular,  the  relationship 
between  GDP  and  emissions  used  to  be  considered  as  linear,  but 
however,  in  certain  cases  this  is  no  longer  true.  The  environmental 
Kuznets  curve  (EKC)  was  firstly  defined  in  the  early  1990s  almost 
simultaneously  by  Shafik  and  Bandyopadhyay  [3],  Panayotou  [11], 
and  Grossman  and  Krueger  [12],  who  stated  that  EKC  refers,  by 
analogy,  to  the  inverted  U-shaped  relationship  between  the  level 
of  economic  development  and  the  degree  of  income  inequality 
[13], 1  According  to  the  EKC  hypothesis,  the  relationship  between 
income  per  capita  and  some  types  of  pollution  is  approximately  an 
inverted  U.  This  behavior  states  that  as  the  GDP  per  capita 
growths,  environmental  damage  increases,  reaches  a  maximum, 
and  then  declines.  The  reason  for  this  behavior  is  that  when  GDP 
reaches  a  certain  threshold,  the  economy  moves  into  a  different 
regime,  where  the  rate  of  emissions  with  respect  to  income  can  be 
reduced  with  respect  to  the  initial  regime.  In  the  initial  stage,  as  in 
the  developing  countries,  C02  emissions  scale  with  the  size  of  the 


1  Exhaustive  surveys  were  made  by  Stern  [14],  Dinda  [15],  and  more  recently 
by  Pasten  and  Figueroa  [16], 


economy  because  the  industries  are  relatively  rudimentary,  unpro¬ 
ductive,  and  polluting.  In  the  second  stage,  the  impact  of  the 
economy  in  environmental  degradation  is  reduced  through  the 
structure  and  composition  effect,  because  the  economy  growth 
induces  structural  changes.  In  particular,  that  happens  as  an 
agricultural  based  economy  shifts  into  a  manufacturing  services 
based  economy.  Finally,  the  third  stage  appears  when  nations 
invest  intensively  in  research  and  development  and  the  dirty  and 
obsolete  technologies  are  replaced  by  clean  ones.  At  this  point  the 
pollution  starts  to  decrease  as  a  function  of  the  GDP.  The  different 
phases  of  the  EKC  are  depicted  schematically  in  Fig.  1. 

The  proposal  of  the  present  work  goes  a  step  further,  and 
intends  to  see  under  which  conditions  a  given  country  could 
approach  the  fulfillment  of  the  EKC  hypothesis  in  the  medium 
term.  To  this  end  we  use  a  model  which  is  based  on  an  extension 
of  the  Kaya  identity2  [17]  and  on  a  GDP  formation  approach  which 
links  renewable  energy  and  income  [18],  The  methodology  was 
first  presented  by  Robalino-Lopez  et  al.  [1,2]  and  was  applied  to 
Ecuador,  here  it  will  be  applied  to  Venezuela. 

One  of  our  aims  is  to  estimate  C02  emissions  in  Venezuela  during 
the  period  2011-2025,  using  a  model  already  presented  in  Robalino- 
Lopez  et  al.  [1,2],  under  four  different  scenarios  which  constrain  GDP, 
productive  sectoral  structure,  energy  intensity,  and  energy  matrix. 
Moreover,  we  will  tiy  to  disentangle  how  these  factors  affect  C02 
emissions  and  to  know  whether  and  when  EKC  hypothesis  could  be 
accomplished  in  Venezuela  in  the  medium  term. 

This  study  could  serve  as  a  guide  for  other  oil-producing 
developing  countries,  because  single  country  studies  can  help 
policy  makers  to  create  comprehensive  policies  to  control  envir¬ 
onmental  degradation. 

The  paper  is  structured  as  follows:  Section  2  presents  an  over¬ 
view  of  the  literature,  Section  3  provides  an  overview  of  the  country, 
Section  4  introduces  the  model  and  the  methodology,  Section  5  is 
devoted  to  discussion  of  empirical  results,  and  finally,  Section  6 
provides  the  summary,  conclusions  and  policy  implications. 


2.  Literature  review 

Nowadays,  the  analysis  of  the  relationship  between  economic 
growth,  energy  consumption  and  C02  emissions  is  a  very  relevant 
topic  in  the  literature  on  sustainable  development  (see  Table  1). 
Pollution  due  to  economic  growth  is  one  of  the  most  important 
empirical  relationships  that  have  been  tested  and  studied  in  the 
ecological  literature  for  many  countries,  both  developing  and 
developed  countries.  The  literature  can  be  divided  into  three  main 
lines  of  study:  first,  the  study  of  the  relationship  GDP-energy 
consumption  or  GDP-energy-C02  emissions,  including  the  study 
of  the  causality  relationships;  second,  the  study  of  the  different 
aspects  of  the  EKC  hypothesis;  and  finally,  the  use  of  scenarios  to 
be  able  to  conduct  forecast  calculations  of  C02  emissions  in  a 
forthcoming  period.  In  many  works  these  lines  of  research  overlap. 
In  Table  1  we  summarize  the  set  of  references  that  will  be 
analyzed  in  detail  below. 


2  C02  emissions  of  a  given  country  are  broken  down  into  the  product  of  four 
factors:  carbon  intensity,  energy  intensity,  economic  rent,  and  population. 
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Regarding  the  first  line  of  study,  in  [19]  the  authors  studied  the 
inequality  in  C02  emissions  across  countries  and  the  relation¬ 
ship  with  the  income  inequality  for  the  period  1971-1999.  They 
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Fig.  1.  Schematic  plot  of  the  U  inverted  relationship  between  GDP  per  capita  and 
C02  emission  per  capita:  stage  1  corresponds  to  a  rapid  growth  of  the  emission, 
stage  2  to  the  stabilization  phase,  and  stage  3  to  the  reduction  of  the  C02  emission 
as  the  income  increases. 


concluded  that  inequality  in  C02  is  mostly  explained  by  the 
inequality  in  the  per  capita  income.  The  relationship  between 
energy  consumption  (electricity  in  this  case)  from  renewable  and 
non-renewable  sources  and  economic  growth  for  18  Latin  Amer¬ 
ican  countries  during  the  period  1980-2010  was  studied  by  Al- 
mulali  et  al.  [20],  The  authors  used  the  Vector  Error-Correction  and 
found  that  renewable  electricity  consumption  is  more  significant 
than  non-renewable  electricity  consumption  in  promoting  eco¬ 
nomic  growth  in  the  long  and  the  short  term.  A  similar  study  was 
carried  out  by  Cowan  et  al.  [21  ]  for  the  BRICS  countries  during  the 
period  1990-2010.  The  authors  explored  the  relationship  between 
electricity  consumption,  economic  growth  and  C02  emissions  and 
they  found  that  the  neutrality  hypothesis  holds  for  Brazil,  India 
and  China,  a  feedback  hypothesis  for  Russia  exits,  a  one  way 
Granger  causality  running  from  GDP  to  C02  in  South  Africa 
appears  and  a  reverse  relationship  from  C02  to  GDP  in  Brazil  is 
found.  In  [22]  the  author  studies  the  C02  emission  flows  and 
energy  use  in  Argentina,  Brazil,  Colombia,  Mexico  and  Venezuela 
from  1971  to  2001.  Results  showed  that  changes  in  C02  intensities 
are  similar  in  all  these  countries,  however  the  energy  use  slightly 
changes  from  country  to  country,  leading  to  small  differences  in  the 


Table  1 

Summary  of  recent  studies  that  analyzed  the  GDP  growth-energy-C02  emissions  relationship. 


Author 

Relationship  Region 

Methodology 

Period 

Outcomes 

Esteve  and  Tamarit 

[27,28] 

C02-GDP  Spain 

EKC 

1857- 

2007 

Existence  of  EKC 

Soytas  and  Sari  [30] 

Energy-GDP  Argentina 

Causality  relationship 

1950- 

1990 

Bi-directional  causality 

Alam  et  al.  [36] 

C02-energy-  India 

GDP 

Dynamic  modeling  and  causal 
relationship 

1960- 

1995 

Bi-directional  causality,  energy-C02  in  the 
long-run  but  neither  C02  nor  energy  causes 
movements  in  income 

Soytas  et  al.  [29] 

Energy-GDP  Turkey 

Cointegration 

1960- 

1995 

Unidirectional  causality,  energy ->  GDP 

Chontanawat  et  al. 
[33] 

Energy-GDP  Central  and  South  America, 
Turkey 

Causality  relationship 

1960- 

2000 

Mixed  evidences 

Halicioglu  [34] 

C02-energy-  Turkey 

GDP 

Causality  relationship 

1960- 

2005 

Strong  connection  GDP-C02 

Kuntsi-Reunanen 

[45] 

C02-energy  Selected  Latin  American 
countries 

C02  emission  flows 

1970- 

2001 

No  significant  changes 

Padilla  et  al.  [19] 

C02-GDP  Groups  of  countries 

Non-parametric  estimations 

1971- 

1999 

GDP  inequality -*C02  inequality 

Shahbaz  et  al.  [39] 

C02-energy-  Pakistan 

GDP 

Cointegration,  Granger  causality  &  EKC 
analysis 

1971- 

2009 

Unidirectional  causal  relationship, 

GDP-C02.  EKC  hypothesis  supported 

Lee  [31] 

Energy-GDP  South  America 

Causality  relationship 

1975- 

2001 

Relationship  energy-* GDP 

Lise  [32] 

C02-energy-  Turkey 

GDP 

Decomposition  analysis 

1980- 

2003 

Decoupling  between  C02  emissions  and  economic 
growth 

Narayan  and 

Narayan  [25] 

C02-GDP  43  developing  countries 

EKC  analysis 

1980- 

2004 

35%  of  the  countries  show  EKC  evidences 

Jaunky  [5] 

C02-energy-  36  high-income  countries 
GDP 

EKC  analysis 

1980- 

2005 

GDP-*C02  EKC  evidence  found 

Ozturk  and  Acaravci 

[35] 

Energy-GDP  Albania,  Bulgaria,  Hungary, 
Romania 

Causality  relationship 

1980- 

2006 

Bi-directional  causality  in  Hungary.  No  causal 
relationships  for  Albania,  Bulgaria  and  Romania 

Shahbaz  et  al.  [40] 

C02-energy-  Romania 

GDP 

Cointegration,  EKC  analysis 

1980- 

2010 

Long  run  relationship  between  economic  growth, 
energy  and  C02.  EKC  evidence  is  found 

Al-mulali  et  al.  [20] 

Energy-GDP  Latin  America 

Cointegration,  Granger  causality 

1980- 

2010 

Long  run  relationship  between  renewable  energy 
and  GDP  growth 

Robalino-Lopez 
et  al.  [1] 

C02-energy-  Ecuador 

GDP 

System  dynamics  modeling  and 
scenario  analysis 

1980- 

2025 

Strong  connection  GDP-C02 

Robalino-Lopez 
et  al.  [2] 

C02-energy-  Ecuador 

GDP 

System  dynamics  modeling  and  EKC 
analysis 

1980- 

2025 

EKC  not  fulfilled 

Cicea  et  al.  [24] 

C02-GDP  European  Union 

Indicator  analysis 

1990- 

2008 

Mixed  evidences 

Cowan  et  al.  [21  ] 

C02-energy-  BRICS  countries 

GDP 

Granger  causality 

1990- 

2010 

Evidences  of  GPD-»C02 

Herrala  et  al.  [23] 

C02-GDP  170  countries 

Stochastic  cost  analysis 

1997- 

2007 

Progress  in  global  efficiency 

Ibrahim  et  al.  [26] 

C02-GDP  69  countries 

GMM  estimators 

2000- 

2008 

Mixed  evidences 

Feng  et  al.  [42] 

Energy,  C02  City  of  Beijing 

System  dynamics  modeling 

2005- 

2030 

Economic  growth  mainly  modulate  C02  emissions 
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evolution  of  C02  emissions.  Hernia  and  Goel  [23]  examine  global  C02 
efficiency  using  a  stochastic  cost  frontier  analysis  of  170  countries 
during  the  period  1997-2007.  The  estimates  show  rather  large 
differences  in  efficiency  levels  and  efficiency  changes  across  countries. 
In  [24]  the  authors  study  the  European  Union  countries  from  1990  to 
2008  using  econometric  models  based  on  the  Kaya  identity.  The 
authors  defined  a  new  environmental  efficiency  index. 

The  study  of  the  EKC  hypothesis  goes  a  step  further  than  the 
above  works,  because,  besides  studying  the  GDP-C02  emissions 
relationship,  it  also  considers  the  EKC  fulfillment.  In  [25]  carbon 
dioxide  emissions  and  economic  growth  have  been  studied  for  a 
panel  of  developing  countries  (43  countries)  in  the  period  1980- 
2004  and  they  found  that  in  the  35%  of  the  sample  the  C02 
emissions  diminished  in  the  long  run,  which  confirm  that  these 
countries  approach  the  sloping  down  part  of  the  EKC.  A  similar 
study  was  conducted  for  a  panel  of  69  countries  by  Ibrahim  and 
Law  [26]  in  the  period  2000-2008  using  the  so-called  generalized 
method  of  moments.  The  paper  finds  evidences  substantiating  the 
presence  of  the  EKC.  Jaunky  [5]  tried  to  test  the  EKC  hypothesis  for 
35  high-income  countries  for  the  period  1980-2005.  He  applied 
several  panel  data  unit  root  and  co-integration  tests.  The  empirical 
tests  provide  evidences  of  an  EKC  for  various  countries.  Esteve  and 
Tamarit  [27,28]  studied  the  EKC  hypothesis  for  the  case  of  Spain 
during  the  period  1857-2007.  The  authors  found  evidences  point¬ 
ing  to  the  existence  of  the  EKC  in  Spain.  A  review  of  the  literature 
shows  a  long-run  relationship  linking  income  and  emissions 
growth  and  recent  research  has  assessed  this  relationship  employ¬ 
ing  cointegration  techniques.  The  empirical  evidence  suggests  that 
GDP  and  pollution  levels  may  be  jointly  determined,  so  that  any 
constraint  put  on  energy  consumption  to  help  the  reduction  of  the 
emissions  will  have  effects  on  economic  growth.  Some  authors 
such  as  Soytas  et  al.  [29],  Soytas  and  Sari  [30],  Lee  [31],  Lise  [32], 
Chontanawat  et  al.  [33],  Halicioglu  [34],  Ozturk  and  Acaravci  [35], 
Alam  et  al.  [36],  Fosten  et  al.  [37],  and  Shahbaz  et  al.  [39,40], 
among  others,  use  cointegration  procedures  to  examine  the 
relationship  between  C02  and  income.  All  these  studies  focused 
on  the  analysis  of  past  evidences. 

In  regard  to  the  third  branch  of  study,  Bautista  [41  ]  conducted 
an  scenario  analysis  for  Venezuela  up  to  2050  paying  special 
attention  to  sustainable  scenarios.  He  showed  that  Venezuela 
has  all  the  resources  needed  to  achieve  a  sustainable  development 
in  the  power  generation  sector  and  that  the  easiest  way  to  reduce 
the  emissions  is  to  improve  the  energy  efficiency.  In  [42]  the 
authors  developed  a  system  dynamics  model  to  describe  the 
evolution  of  energy  consumption  and  C02  emissions  in  the  city 
of  Beijing  from  2005  to  2030.  They  conducted  a  sensitivity  analysis 
showing  that  the  economic  growth  and  the  increase  of  the 
population  are  the  two  main  variables  that  modulate  the  value 
of  C02  emissions.  In  [1]  the  authors  developed  a  model  based  on 
the  Kaya  identity  and  on  the  formation  of  the  GDP  that  depends  on 
the  use  of  renewable  energies.  The  model  was  applied  to  Ecuador 
for  the  period  1980-2025.  The  authors  carried  out  a  scenario 
analysis  and  concluded  that  in  order  to  control  the  C02  emission 
under  scenarios  of  rapid  growth  of  the  GDP  there  are  two 
complementary  strategies,  to  improve  the  energy  efficiency  and/ 
or  to  promote  the  use  of  renewable  energy. 

It  deserves  to  be  mentioned  in  Ref.  [2],  because  the  authors 
study  the  EKC  fulfillment  not  with  past  data,  but  in  a  forthcoming 
period.  They  studied  Ecuador  from  1980  to  2025,  starting  with  the 
model  proposed  in  [1 J.  The  main  conclusion  of  this  work  is  that 
Ecuador  will  enter  into  the  second  stage  of  the  EKC  in  the  coming 
years  if  the  country  follows  one  of  the  scenarios  where  the  energy 
efficiency  and  the  renewable  energies  are  promoted. 

After  this  review,  it  is  clear  the  lack  of  studies  in  the  literature 
on  the  EKC  hypothesis  in  a  forthcoming  period,  in  South  America, 
in  general,  and  in  Venezuela,  in  particular. 


3.  Overview  of  the  study  area 

Venezuela  has  abundant  energy  resources  that  have  played  a 
key  role  in  defining  the  nation's  economic  development  road. 
The  extraction  of  oil  has  formed  the  backbone  of  the  Venezuelan 
economy  since  1920s.  The  economy  of  Venezuela  is  largely  based 
on  petroleum  and  manufacturing  sectors.  In  2010  revenue  from 
petroleum  exports  accounted  for  approximately  18%  of  the  coun¬ 
try's  GDP  and  roughly  93%  of  total  exports  [43],  Hence,  the 
Venezuelan  GDP  is  highly  correlated  with  the  price  of  the  oil, 
therefore  it  suffers  large  fluctuations;  during  the  collapse  of  oil 
prices  in  the  1980s,  the  economy  contracted,  however  it  grew  an 
average  of  10.5%  between  2004  and  2008,  but  shrank  an  average  of 
2.3%  from  2009  to  2010.  The  country  also  has  one  of  the  highest 
inflation  rates  in  the  world,  averaging  29.1%  in  2010  [43],  Accord¬ 
ing  to  Central  Bank  of  Venezuela  (BCV-Banco  Central  de  Vene¬ 
zuela),  the  government  received  around  325  billion  USD  through 
oil  production  and  exports  from  1998  to  2008  [44], 

Venezuela  had  the  highest  energy  consumption  per  capita  and 
the  highest  energy  intensity  in  Latin  America.  According  to  IEA  the 
primary  energy  use  is  about  20%  higher  than  it  would  be  without 
subsidies  [45].  Venezuela  has  an  oil  production  quota  that  cannot 
be  exceeded  because  it  is  an  active  member  of  OPEC;  therefore, 
a  high  domestic  consumption  of  oil  penalizes  the  exports  and  the 
country  revenues. 

The  other  side  of  the  high  oil  consumption  is  the  high  level  of 
C02  emissions  of  the  country;  in  2009  Venezuela  emitted 
6.5  metric  tons  of  C02  per  capita.  This  makes  Venezuela  the  largest 
C02  emitter  per  capita  in  Latin  America,  very  far  though  from 
Qatar,  the  world's  largest  C02  emitter  with  44  metric  tons  per 
capita.  It  is  expected  that  the  present  social  and  economic 
development  in  the  region  could  significantly  increase  the  Venezuelan 
level  of  emissions. 

Since  2000,  the  annual  average  primary  energy  production  of 
Latin  America  and  the  Caribbean  reached  447  million  tonnes  of  oil 
equivalent  (Mtoe).  This  production  was  made  of  70%  of  carbon 
fossil  fuels  (46%  crude  oil,  20%  gas  and  4%  coal),  the  rest  divided 
between  biomass  and  garbage  (18%),  hydroelectricity  (10%), 
nuclear  (1%),  and  geothermic,  solar  and  other  renewable  sources 
of  energy  (1%)  [43],  The  energy  mix  in  Venezuela  is  different  for 
this  pattern,  with  a  higher  relative  importance  of  carbon  fossil 
fuels,  which  reached  89%  of  the  total  primary  energy  production 
(47%  of  oil  and  42%  of  gas),  with  10%  of  hydroelectricity  and  1%  of 
biomass  and  garbage  by  2010  [43,44], 

The  Venezuelan  government  has  identified  potential  sources  of 
renewable  energy  that  would  enable  a  reduction  of  the  consumption 
of  fossil  fuel  in  the  near  future  [38,41],  At  present,  Venezuela  has  a 
hidro  power  complex  at  Carol  river  with  and  annual  energy  produc¬ 
tion  of  130  TWh,  but  the  potential  energy  yield  of  this  complex  could 
reach  246  TWh.  Regarding  wind  and  solar  power,  nowadays  the 
installed  power  is  negligible,  but  the  potential  resources  are  huge: 
1038  TWh/year  and  1465  TWh/year  of  wind  and  solar  power,  respec¬ 
tively.  In  summary,  Venezuela  has  available  renewable  energetic 
resources  that  can  significantly  decrease  the  contribution  of  fossil 
energies  to  the  energy  matrix  of  the  country  [46], 


4.  Model  and  methodology 

4.2.  Formulation  of  the  model 

The  model  to  calculate  C02  emissions  from  fossil  energy 
corresponds  to  a  nexus  relationship,  which  is  an  extension  of  the 
original  Kaya  identity  where  we  disaggregate  in  type  of  fuel  and 
kind  of  industrial  sector  (see  Refs.  [1,2]  for  a  complete  description 
of  the  model).  The  amount  of  C02  emissions  from  industry  and 
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other  energy  uses  may  be  studied  quantifying  the  contributions  of 
five  different  factors:  scale  of  the  economy,  industry  activity  mix, 
sectoral  energy  intensity,  sectoral  energy  mix,  and  C02  emission 
factors.  Moreover,  we  consider  different  sub-categories  concerning 
industrial  sectors  and  fuel  type.  The  C02  emissions  can  be  written 
as 

c  =  zc,j  =  §■  f  % =  QZ Si  ■  Eh  ■  My  ■  Ulp  (1 ) 

where  C  is  the  total  C02  emissions  (in  a  given  year);  Cg  is  the  C02 
emission  arising  from  fuel  type  j  in  the  productive  sector  i  (note 
that  the  index  i  runs  over  four  productive  sectors  and  the  index  j 
over  seven  type  of  energy  sources);  Q  is  the  total  GDP  of  the 
country;  Q,  is  the  GDP  generated  by  the  productive  sector  i;  E,  is 
the  energy  consumption  in  the  industrial  sector  i;  Eg  is  the 
consumption  of  fuel  j  in  the  industrial  sector  i;  S,  (Q,/Q)  is  the 
share  of  sector  i  to  the  total  GDP;  the  energy  intensity  of  sector  i  is 
given  by  £/,■  (E ,/Q,);  the  energy  matrix  is  given  by  Mg  (Ey/Ef)  and 
the  C02  emission  factor  by  Ug  (Cy/Ey).3  The  definition  of  the 
productive  sector  division  and  the  classification  of  type  of  fuels 
comply  the  standard  separation  of  the  International  Standard 
Industrial  Classification  [47]  and  the  International  Energy  Agency 

[48] ,  respectively. 

It  is  also  of  interest  to  write  up  how  to  calculate  the  total 
energy  in  terms  of  the  GDP, 

E=ZEy  =  QZSiEI,My  (2) 

ij  ‘j 

and  the  expended  energy  of  eveiy  kind  of  fuel, 

Ej  =  QZS,  ■  Eh  ■  My.  (3) 

The  subsequent  data  analysis  and  the  pre-processing  of  the 
time  series  were  performed  using  the  Hodrick-Prescott  (HP)  filter 

[49] ,  which  allows  isolation  of  outliers  (economic  crises,  random 
behavior  of  markets,  etc.)  of  the  time  series  under  study.  After  that, 
it  is  possible  to  get  the  trend  component  of  a  time  series  and  to 
perform  estimations  that  are  more  appropriate.  The  smoothing 
parameter,  X,  of  the  filter,  which  penalizes  acceleration  in  the  trend 
relative  to  cycle  component,  is  taken  equal  to  100.  Most  of  the 
business  cycle  literature  uses  this  value  for  the  X  parameter,  as  has 
been  suggested  by  Hodrick  and  Prescott  [49], 

4.2.  Economic  submodel 

This  work  presents  as  a  key  point  the  explicit  inclusion  of  the 
effect  of  renewable  energy  on  the  GDP,  allowing  to  establish  a  link 
between  income  and  C02  emissions  [50],  In  this  work  we  consider 
that  renewable  energy  can  increase  GDP  through  import  substitu¬ 
tion  of  energy  which  has  direct  and  indirect  effects  on  increasing 
GDP  and  trade  balance. 

Closely  following  Chien  and  Hu  [18]  we  use  the  expenditure 
approach  to  form  the  GDP, 

( l  =  Ca  +  I  +  G+TB ,  (4) 

where  Ca  is  the  final  household  consumption  expenditure,  /  the 
gross  domestic  capital  formation,  G  is  the  general  government 
final  consumption  expenditure,  and  TB  the  trade  balance  (exports 
minus  imports). 

On  the  other  hand,  according  to  Chien  and  Hu  [18],  the  variable 
G  should  be  removed  from  the  model  estimation  to  avoid  multi- 
collinearity.  The  system  of  theoretical  GDP  formation  model  is 


3  Throughout  this  paper,  as  a  convention,  we  will  always  refer  to  the  industrial 
sector  with  the  i  index  and  to  the  type  of  energy  source  with  the  j  index. 


composed  by  the  following  equations: 

Q_=  Ci\  ■  I+a2  ■  TB  -j-  tq  •  Ca  4-  Q4  ■  Ejmp  +  05-  RN  +  £j  (5) 

I  =  bi  ■  RN+b2  ■  Ca+e2  (6) 

TB  =  Ci  •  E[mp  4-  c2  ■  RN 4- €3  (7) 

Eimp  =  dl  ■  RN+£4  (8) 

Ca  =f  1  ■  Ejmp+f2  ■  TB  +  £5  (9) 


where  Eimp  is  the  energy  import,  RN  is  the  renewable  energy  and 
ei...e5  are  residuals.  The  coefficients  at,  b„  q,  d„  and/)  appearing  in 
these  equations  are  determined  using  the  Seemingly  Unrelated 
Regression  (SUR)  with  the  smoothed  dataset  of  the  period  1980- 
2010.  The  data  used  to  calibrate  the  model  corresponds  to  the 
period  1980-2010  and  was  extracted  from  the  official  dataset  of 
the  country.4  Table  2  shows  the  results  after  fitting  the  smoothed 
series  of  data. 

In  Eq.  (5),  GDP  is  positively  influenced  by  capital  formation, 
trade  balance,  and  consumption  (cq,  a2,  and  a3  are  positive,  as  seen 
in  Table  2).  Chien  and  Hu  [51]  suggested  that  energy  inputs  may 
increase  GDP,  as  such,  energy  imports  and  renewable  energies  are 
included  as  well  in  Eq.  (5).  In  the  case  of  Venezuela,  our  results 
confirm  that  the  import  of  energy,  Eimp,  has  a  small,  but  positive 
influence  (a4>0)  on  the  GDP,  which  is  somehow  surprising:  a 
larger  energy  import  implies  a  larger  income.  Finally,  renewable 
energies  also  have  a  positive  effect  on  the  GDP  (a5  >  0),  confirming 
the  theory  presented  in  Chien  and  Hu  [51],  This  is  not  a  general 
property  and  strongly  depends  on  the  region  under  study.  In 
particular,  in  the  case  of  Ecuador  a5  <0  (see  [1.2]). 

In  Eq.  (6)  capital  formation  is  influenced  by  renewable  energy, 
thus  increasing  the  use  of  renewable  energies  will  result  in 
business  expansion  and  thus  capital  could  be  accumulated  in  long 
term  (fq  and  b2  are  positive  in  Table  2).  In  Eq.  (7)  a  reduction  of 
energy  imports  and  an  increase  of  the  use  of  renewable  energy  has 
a  positive  influence  in  trade  balance  (ci  is  negative  while  c2 
positive  in  Table  2;  note  that  Eimp  <  0  for  net  energy  exports,  as 
it  is  the  case  of  Venezuela).  The  theory  proposed  by  Domac  et  al. 

[50]  suggests  that  the  use  of  renewable  energy  results  in  import 
substitution  and  thus  trade  balance  will  increase  with  the  use  of 
renewable  energies.  This  is  particularly  true  in  the  case  of 
Venezuela  where  oil  extraction  is  limited  due  to  the  OPEC  quota. 
Although  Venezuela  is  a  net  exporter  of  fossil  energy,  their  fossil 
energy  consumption  could  be  supplied  by  renewable  energy,  thus, 
diversifying  its  energy  matrix  and  reducing  the  emissions.  Hence, 
energy  imports  will  be  reduced  by  increasing  the  use  of  renewable 
energy  (in  Eq.  (8)  coefficient  dj  is  negative  which  confirms  this 
assumption).  In  Eq.  (9),  according  to  international  trade  theories, 
the  domestic  price  of  goods  increases  as  the  same  kind  of  goods 
are  exported,  while  it  decreases  as  the  same  kind  of  goods  are 
imported.  Thus,  trade  balance  (/1  >  0)  influences  consumption 
through  changes  in  domestic  prices,  as  well  as  imports  of  energy 
( f2  <  0)  influence  domestic  energy  prices  and  therefore  increase 
the  household  consumption.  As  a  result,  consumption  of  energy- 
related  products  is  also  affected. 

4.3.  Energy  consumption  and  productive  sectoral  structure  submodel 

Energy  consumption  refers  to  the  use  of  primary  energy  before 
transformation  into  any  other  end-use  energy,  which  is  equal  to 


4  Data  is  from  World  Bank  [43],  and  International  Energy  Agency  [48], 
Economic  official  data  set  used  is  given  in  constant  2005  PPP  (purchasing  power 
parity)  international  dollars  [43],  In  the  rest  of  this  paper  GDP-PPP  will  be  referred 
to  only  as  GDP  and  2005  PPP  international  dollars  as  USD,  for  brevity. 
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Table  2 

Estimated  coefficients  for  the  GDP  formation  equations  using  the  SUR  technique 
(see  Eqs.  (5)-(9)). 


Variable 

Of 

/ 

TB 

F 

•-imp 

c 

(a,) 

( bi ) 

(Q) 

(d,) 

(fi) 

/a 

1  4g*** 

TBa 

(152.06) 

0.88*** 

1.39*** 

<? 

(22.66) 

0.85*** 

3.29  io-2*** 

(4.96) 

Eimpb 

(61.39) 

6.16  10~5*** 

(0.39) 

-4.36  10-6 

-1.48  10~4 

RNb 

(14.40) 

2.19  10~3  *** 

6.45  io-3*** 

(-0.17) 

4.04  io-3*** 

-11.95*** 

(-1.41) 

(8.01) 

(7.30) 

(12.30) 

(-8.83) 

a  a  I.  TB,  C  given  in  109  USD. 
b  Eimp  and  RN  given  in  103  ktoe. 

***  Significance  at  the  1%  level  and  numbers  in  the  parentheses  are  t-statistics. 


the  local  production  of  energy  plus  imports  and  stock  changes, 
minus  the  exports  and  the  amount  of  fuel  supplied  to  ships  and 
aircrafts  engaged  in  international  transport.  It  is  given  in  ktonnes 
of  oil  equivalent  (ktoe).  Energy  intensity  is  defined  as  the  ratio  of 
energy  consumption  and  GDP  [43]. 

In  this  work  we  consider  four  sectors  to  define  the  productive 
sectoral  structure  which  according  to  the  International  Standard 
Industrial  Classification  of  all  economic  activities  [47]  are:  (1)  pri¬ 
mary  sector  (Sector  1:  agriculture,  forestry,  and  fishing),  (2)  sec¬ 
ondary  sector-A  (Sector  2:  manufacturing),  (3)  secondary  sector-B 
(Sector  3:  mining,  quarrying,  electricity,  gas,  water  supply,  and 
construction),  and  (4)  tertiary  sector  (Sector  4:  services,  trade, 
residential,  and  transportation).  These  are  represented  inside  the 
model  by  their  contribution  to  the  country’s  economy  (S,),  by  their 
energy  intensity  (£/,),  and  by  their  energy-fuel  mix  (My).  Index  i 
runs  over  each  sector  of  the  productive  sectoral  structure,  while 
index  j  runs  over  seven  type  of  fuels:  (1)  natural  gas,  (2)  liquefied 
petroleum  gases  (LPG),  (3)  motor  gasoline,  (4)  gas/diesel  oil, 
(5)  fuel  oil,  (6)  petroleum  coke,  and  (7)  renewable,  alternative, 
and  nuclear  (not  present  in  Venezuela)  energy. 

4.4.  C02  intensity  and  energy  matrix  submodel 

C02  intensity  (C02mt)  of  a  given  country  corresponds  to  the 
ratio  of  C02  emissions  and  the  total  consumed  energy  written  in 
terms  of  mass  of  oil  equivalent  (C02i„r  =  HijQj/EijEy)-  The  value  of 
the  CO 2int  in  a  given  year  depends  on  the  particular  energy  mix 
during  that  year.  My  gives  the  energy  matrix,  but  it  is  more 
convenient  to  sum  over  the  different  sectors  and  aggregate  the 
fossil  fuel  contributions,  therefore,  we  define  Mj  =  £iEy/£ijEy  for 
each  type  of  fuel  (j=  1,...,7).  Hence,  is  the  share  of 

expended  fossil  energy,  while  M7  is  the  share  of  used  renewable 
energy  in  the  country. 

In  order  to  simplify  the  description,  in  this  paper  we  assume 
that  M7  does  not  contribute  to  the  C02  emissions.  Following  the 
methodology  recommended  by  the  IPCC,  that  is,  the  Reference 
method,  the  approach  of  the  first  level  (Tier  I)  for  the  fossil  energy 
mix  was  used.5 

4.5.  Model  equations  and  causal  diagram 

In  the  equations  presented  along  this  section  we  have  seen  how 
the  model  is  split  in  two  different  parts:  energy  and  productive 


5  The  emission  factors,  ! are  taken  from  the  IPCC  methodology  to  estimate 

C02  emission  of  each  fuel  [8], 


sectoral  submodel,  Eqs.  (2)  and  (3),  and  economic  submodel, 
Eqs.  (5)-(9)).  In  the  first  case,  the  energy  and,  in  particular,  the 
amount  of  renewable  energy,  Ej  =  1,  are  calculated.  In  the  second 
one,  the  value  of  the  GDP  is  calculated  in  terms  of  its  components, 
one  of  which  is  the  renewable  energy.  These  two  parts  are  coupled 
through  the  renewable  energy  term  which  generates  a  feedback 
mechanism  that  in  the  case  of  Venezuela  is  positive  (a5  >  0,  see 
Table  2).  In  Fig.  2  we  present  the  schematic  view  of  the  model, 
where  the  feedback  mechanism  is  highlighted.  This  way  of  pre¬ 
senting  the  model  is  extremely  useful  because  it  allows  to  see  the 
driving  forces  of  C02  emissions  in  a  hierarchical  way,  showing  the 
causality  relationship  between  the  different  variables.  It  can  be 
observed  that  the  C02  emitted  into  the  atmosphere  has  several 
connections  with  the  variables  of  the  model:  economic  growth, 
productive  sectoral  structure,  energy  consumption,  and  energy 
matrix. 

A  more  quantitative  way  of  presenting  how  the  different 
variables  are  extrapolated  is  through  the  difference  equations  that 


should  be  solved: 

Q(  f)  =  a,  /(t)  +  a2  TB(  t )  +  a3  Ca(  t)  +  a4Eimp(t )  +  a5RN(  t  - 1 ),  (10) 

Ey(t)  =  £Sft)  ■  Elft)  ■  Mij(t )  ■  CDP(t),  (1 1 ) 

i 

RN(t)  =  E7(t),  (12) 

y(t)  =  y(t - 1)  ■  (l  +  ry),  (13) 


where  Sft),  EI,{t),  My(t),  I(t),  TB(t ),  Ca(t),  and  Eimp(t)  evolve  follow¬ 
ing  Eq.  (13)  while  the  parameters  a,  have  constant  values.  Note 
that  index  j  runs  over  the  type  of  energy  sources,  while  i  on  the 
industrial  sectors;  j=  7  corresponds  to  renewable  and  alternative 
energy.  t=0  corresponds  to  the  reference  year  2010  and  f  is  given 
in  number  of  years  since  2010.  The  value  of  ry  is  fixed  through  the 
definition  of  the  used  scenario.  In  the  case  of  the  BS  scenario  one 
should  use  a  value  of  ry  that  depends  on  the  time  [52,53].  In  this 
case  ry  roughly  corresponds  to  the  yearly  average  increase  over  the 
period  1980-2010. 

The  feedback  mechanism  is  provided  through  the  inclusion  of 
KN(f-l)  in  the  calculation  of  the  GDP  (10).  As  a5  >  0  (see  Table  2) 
the  feedback  mechanism  is  positive.  This  fact  induces  an  extra 
increase  of  the  GDP  for  the  SC-3  and  SC-4  scenarios  with  respect  to 
SC-2  (see  Section  4.7)  for  increasing  the  use  of  renewable  energy. 
In  general,  any  increase  of  the  term  £,• S,-(t)  •  Elft)  ■  My(t)  for  j=T 
will  induce  an  increasing,  though  moderate,  of  the  GDP. 

4.6.  Model  validation  and  verification 

The  official  dataset  from  1980  to  2010  and  the  output  of  the 
model  can  be  compared  to  test  its  robustness  and  reliability.  This 
analysis  can  be  carried  out  calculating  the  mean  absolute  percen¬ 
tage  error  (MAPE),  which  is  defined  as 

MAPE(%)  =  1  £ 

i  =  1 

where,  Data „  Mode/,,  and  N  are  the  real  data,  the  calculated  values, 
and  the  number  of  data,  respectively.  MAPE  is  commonly  used  to 
evaluate  cross-sectional  forecasts  [54], 

MAPE  values  for  the  main  variables  of  the  model  are  shown  in 
Table  3.  The  values  are  rather  modest  and  therefore  support  the 
robustness  of  the  model.  To  understand  the  differences  between 
the  real  data  and  the  output  of  the  model  for  the  historical  period, 
one  has  to  take  into  account  that  prior  to  use  the  data  they  are  first 
smoothed  through  the  HP  filter.  On  the  other  hand,  the  description 
of  the  GDP  implies  a  linear  regression  on  smoothed  time  series. 
Finally,  to  calculate  the  C02  emissions  we  use  generic  emission 
factors,  which  generates  additional  discrepancies  when  comparing 


Dataj—Modelj 


Data,- 


x  100, 


(14) 
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Fig.  2.  Causal  diagram  of  the  model.  Continuous  lines  stand  for  the  relationship  between  variables,  while  dashed  ones  correspond  to  control  terms  (S:  productive  sectoral 
structure,  M:  energy  matrix,  U:  emission  factors).  Bold  line  represents  a  feedback  mechanism. 


Table  3 

MAPE  values  for  the  main  variables  of  the  model. 


Variable 

MAPE  (%) 

GDP 

4.6 

Ffossil 

4.4 

RN 

4.4 

Ftotal 

4.1 

co2 

4.2 

with  real  data  (C02  inventory)  that  were  computed  using  a  fully 
disaggregated  sum  and  specific  emission  factors  for  the  country. 
Finally,  note  that  in  some  cases  part  of  the  discrepancies  are 
canceled  out,  in  particular,  the  MAPE  value  for  C02  is  smaller  than 
for  GDP  and  £fossil,  but  however  C02  is  calculated  using  both 
variables  among  other  ingredients. 


4.7.  Scenarios 

In  this  work  we  will  calculate  the  value  of  C02  emissions  in  the 
medium  term.  As  we  have  seen  in  the  previous  sections,  C02 
emissions  depend  on  several  variables,  therefore  it  is  compulsory 
to  define  scenarios  which  determine  how  GDP,  sectoral  structure, 
energy  intensity,  and  energy  mix  will  evolve  in  the  coming  years. 

To  define  the  scenarios  we  will  use  as  a  guide  next  goals:  Goall, 
by  2025  GDP  will  be  substantially  higher  than  in  2010  through  a 
process  of  industrialization  and  improvement  of  the  productive 
sectoral  structure  of  the  country;  Goal 2,  in  regard  to  Goal 1,  the  use 
of  renewable  energies  will  be  twice  the  one  in  2010;  Goal3,  in 
regard  to  Goal  1  and  Goal2,  energy  efficiency  will  increase  with  a 
reduction  of  the  energy  intensity  (2%  per  year  in  each  productive 
sector)  and  by  changes  in  the  productive  sectoral  structure. 

Taking  into  account  the  latter  goals,  we  propose  four  scenarios 
concerning  the  growth  of  GDP,  the  evolution  of  the  energy  matrix, 
of  the  productive  sectoral  structure,  and  the  improvement  of 
energy  efficiency  for  the  period  2011-2025. 

1.  Baseline  scenario  ( BS ):  GDP,  energy  matrix  and  productive 
sectoral  structure  will  evolve  through  the  smooth  trend  of  the 


period  1980-2010,  extrapolated  to  2011-2025  using  the  geo¬ 
metric  growth  rate  method. 

2.  Doubling  of  the  income  (SC-2  scenario):  GDP  in  2025  will  be 
more  than  double  that  of  2010.  To  generate  this  scenario  a 
constant  annual  growth  of  the  GDP  formation  components 
(/,  TB,  C,  Eimp )  of  4.7%  per  year  between  2011  and  2025  is 
assumed  and  a  structural  change  in  the  productive  sectoral 
structure  will  be  implemented  through  a  stabilization  of  the 
share  of  primary  sector  (Sector  1)  at  5%,  a  consolidation  of  the 
country's  industrial  sector  (Sectors  2  and  3)  increasing  its 
share  from  52%  to  55%  and  a  stabilization  of  the  service  sector 
(Sector  4)  at  40%.  The  rest  of  the  variables  will  evolve  as  in  the 
BS  scenario.  This  scenario  clearly  corresponds  to  a  situation 
where  the  economy  is  growing  rapidly  and  no  mitigation 
measurements  to  reduce  C02  emissions  are  carried  out. 

3.  Doubling  of  the  income  and  of  the  share  of  renewable  energies 

(SC-3  scenario):  The  doubling  of  the  GDP  and  the  change  of  the 
productive  sectoral  structure  as  in  the  SC-2  scenario  are 
considered,  however  the  share  of  fossil  energy,  ,  Mj,  will 

be  reduced  approximately  one  point  per  year,  passing  from  88% 
in  2010  to  around  76%  in  2025  due  to  a  constant  annual  growth 
of  the  renewable  and  alternative  energy  share.  This  scenario 
shows  a  first  measure  of  the  environmental  responsibility  in 
order  to  try  to  reduce  dependence  on  fossil  energy. 

4.  Doubling  of  income,  doubling  of  renewable  energy  share  and 
improvement  in  the  efficiency  of  energy  use  (SC-4  scenario):  The 
doubling  of  the  GDP,  the  change  in  the  productive  sectoral 
structure  and  the  change  of  the  share  of  the  fossil  energy  are 
the  same  as  in  the  SC-3  scenario.  Moreover,  an  improvement  in 
the  efficiency  of  energy  use  is  implemented  with  a  reduction  in 
the  sectoral  energy  intensity  of  1%  in  each  productive  sector, 
relative  to  the  BS  scenario,  which  in  practice  supposes  a  2% 
yearly  reduction.  This  scenario  improves  the  country’s  environ¬ 
mental  responsibility  and  sustainable  development  by  support¬ 
ing  energetic  saving  measures  and  energy  efficiency. 


It  is  important  to  note  that  the  aim  of  this  paper  is  not  to 
perform  a  rigorous  forecast  of  the  GDP,  energy  consumption  or 
energy  intensity,  but,  we  try  to  establish  a  baseline  and  other 
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reasonable  scenarios  that  could  be  useful  as  reference  points  for 
policy  makers  or  further  studies. 

4.8.  EKC  hypothesis  analysis 

To  explore  in  a  quantitative  way  the  relationship  between  C02 
and  GDP,  we  follow  Jaunky  [5]  specification  for  testing  the  EKC 
hypothesis  in  Venezuela.  We  consider  the  next  simplified  equation 
that  relates  C02  and  GDP, 

LC02t  =  p0+p^LCDPt+et,  (15) 

where,  LC02  is  the  natural  logarithm  of  the  C02,  while  LGDP  is  the 
natural  logarithm  of  GDP.  p0  is  a  constant  term,  estimates  the 
C02-GDP  elasticity,  and  e  is  an  error  term.  The  use  of  the  natural 
logarithm  is  needed  to  define  in  a  proper  way  the  elasticity  in 
Eq.  (15).  The  region  of  little  environmental  respect  corresponds  to 
p1>  1  and  therefore  to  the  stage  1  of  Fig.  1.  In  this  region  a  change 
in  GDP  induces  a  large  change  in  C02  emissions,  i.e.,  there  is  a  high 
responsiveness.  In  general,  countries  stay  in  this  phase  in  the  early 
stage  of  its  industrial  development.  If  0  <p.x  <  1,  then  an  increase 
in  the  income  generates  an  increase  less  than  proportionate  in  the 
emissions.  The  country  enters  in  the  second  stage  of  the  EKC  with 
an  environmental  stabilization.  Finally,  <  0  implies  that  any 
increase  in  the  GDP  brings  a  decrease  in  the  C02  emissions.  This  is 
the  final  stage  of  the  EKC  and  the  country  is  moving  into  a 
situation  of  intensive  use  of  sustainable  technology  and  GDP 
increase  is  accompanied  by  reductions  of  C02  emissions. 

5.  Empirical  results  and  discussion 

Below,  the  empirical  outcome  of  this  study  is  presented.  In 
particular,  we  provide  estimations  of  GDP  and  C02  emissions  till 
2025,  as  well  as  the  study  of  the  EKC  hypothesis.  For  completeness, 
some  other  projections  are  also  shown.  Besides  we  carry  out  a 
sensitivity  analysis  where  we  disentangle  the  factors  that  drive 
C02  emissions. 

5.1.  GDP 

We  depict  in  Fig.  3  the  value  of  the  GDP  for  the  four  considered 
scenarios.  These  values  do  not  intend  to  be  a  detailed  projection  of 


Fig.  3.  GDP  of  Venezuela  for  the  period  1980-2025.  Dots  correspond  to  the  official 
dataset. 


the  GDP  but  rather  are  a  consequence  of  the  definition  of  the 
scenarios.  The  starting  point  of  our  projection  is  the  data  of  the 
GDP  for  2010  which  is  316  billion  USD.  In  the  case  of  the  BS 
scenario  the  value  of  the  projected  GDP  for  2025  is  484  billion 
USD,  while  the  SC-2  scenario  reaches  643  billions  USD,  which  is 
more  than  double  that  of  2010;  SC-3  and  SC-4  reach  673  and  672 
billions  USD,  respectively.  In  the  case  of  the  BS  scenario  the  yearly 
growth  of  the  GDP  corresponds  to  the  average  of  the  period  1980- 
2010.  The  other  three  scenarios  have  a  annual  growth  of  the  GDP 
components  of  about  5%  and  the  GDP  in  2025  is  roughly  the 
double  that  in  2010.  The  differences  in  the  value  of  the  projected 
GDP  between  the  three  latest  scenarios  come  from  the  presence  of 
a  positive  feedback  mechanism,  which  creates  an  additional 
growth  of  the  GDP  of  about  4.7%  in  2025  with  respect  to  the 
SC-2  scenario. 


5.2.  Energy 

In  this  section  we  present  the  evolution  of  the  energy  con¬ 
sumption,  disaggregating  in  fossil  and  renewable  energy,  for  the 
four  proposed  scenarios.  In  Fig.  4  we  depict  the  evolution  of  the 
energy  consumption  in  the  upper  panel,  the  projection  of  the  fossil 
energy  consumption  in  the  middle  and  the  renewable  energy  in 
the  lower  panel. 

The  total  energy  consumption  in  2025  in  the  BS  scenario  is 
around  102  Mtoe,  while  130  Mtoe  in  the  SC-2, 135  Mtoe  in  the  SC-3 
and  113  Mtoe  in  the  SC-4  scenarios.  Regarding  the  fossil  energy 
consumption  in  2025,  it  corresponds  to  90  Mtoe,  115  Mtoe, 
105  Mtoe,  and  88  Mtoe  in  the  BS,  SC-2,  SC-3,  and  SC-4  scenarios, 
respectively.  Therefore,  renewable  energy  production  is  12  Mtoe, 
15  Mtoe,  30  Mtoe,  and  25  Mtoe  in  2025  for  the  BS,  SC-2,  SC-3,  and 
SC-4  scenarios,  respectively.  The  contribution  from  renewable 
energy  will  allow  to  increase  the  oil  exports  in  about  10-17  Mtoe 
when  comparing  the  SC-3  and  SC-4  with  the  SC-2  scenario.  The 
bigger  renewable  energy  production  in  the  SC-2  scenario  in 
comparison  with  the  BS  scenario  is  because  the  larger  energy 
consumption  due  to  larger  economic  growth  in  the  SC-2  scenario. 


year 


Fig.  4.  Energy  consumption  of  Venezuela  for  the  period  1980-2025.  Top  panel: 
total  energy  consumption.  Mid  panel:  fossil  energy  consumption.  Bottom  panel: 
renewable  energy  consumption.  Dots  correspond  to  the  official  dataset. 
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5.3.  C02  emissions 

In  Fig.  5  we  depict  the  value  of  C02  emissions  during  the  period 
2011-2025,  under  the  four  considered  scenarios.  In  2025  the 
highest  emission  corresponds  to  the  SC-2  scenario  while  the 
lowest  to  the  SC-4  scenario,  though  very  close  to  the  BS  scenario, 
being  the  SC-3  scenario  in  between  these  two  extreme  situations. 
The  SC-2  scenario  does  not  implement  any  attenuation  measure 
regarding  C02  emission  and  moreover  a  rapid  GDP  increase  is 
considered,  therefore  the  larger  value  of  the  emissions  is  expected, 
reaching  in  2025  a  value  of  306  Mtonnes,  while  the  reference 
scenario  only  reaches  240  Mtonnes.  In  the  SC-3  scenario  the  share 
of  renewable  energy  in  the  energy  mix  is  increased  and  the 
emissions  are  reduced  down  to  280  Mtonnes.  Finally,  for  the  SC- 
4  scenario  besides  the  promotion  of  renewable  energies  a  reduc¬ 
tion  of  the  energy  intensity  is  implemented  reducing,  in  this  case, 
the  emissions  down  to  234  Mtonnes.  Therefore,  there  are  two 
different  strategies  to  reduce  the  emissions,  either  promote 
renewable  energies,  which  achieves  a  reduction  of  26  Mtonnes, 
or  to  improve  the  energy  intensity,  which  accomplish  a  reduction 
of  46  Mtonnes.  The  combination  of  both  strategies  leads  to  a 
reduction  of  72  Mtonnes  with  respect  to  the  more  polluted 
scenario,  which  corresponds  to  a  reduction  of  24%  in  2025. 

5.4.  Sensitivity  analysis 

In  this  section  we  will  present  a  sensitivity  analysis  based 
on  the  logarithmic  mean  Divisia  index  (LMDI)  [55],  This  analysis 
will  allow  us  to  determine  the  relative  importance  of  each  term 
appearing  in  the  C02  emission  formula  (1).  In  fact,  using  this 
method  it  is  possible  to  write  down  the  increase  of  C02  emissions 
relative  to  the  value  of  a  given  year,  and  to  decompose  it  as  the 
product  of  the  factors  corresponding  to  the  different  parts  that 
make  up  the  value  of  C02  emission.  Therefore,  we  can  write  [55] 

Dfot  =  Dact  x  Dstr  X  Dint  x  Dmjx  X  Demj  (16) 

where  Dto[  is  the  value  of  the  C02  emission  relative  to  a  reference 
year,  2010  in  our  case,  Dact  accounts  for  the  increase  of  the  GDP 
relative  to  the  reference  year,  Dstr  is  the  change  of  the  structure 
term  (the  share  of  the  different  sectors  to  the  GDP),  Din[  stands  for 


Fig.  5.  C02  emissions  of  Venezuela  for  the  period  1980-2025.  Dots  correspond  to 
the  official  dataset. 


the  change  of  the  intensity  term  that  refers  to  energy  consump¬ 
tion,  and  Dmix  accounts  for  the  change  in  the  energy  mixing  term.6 

This  analysis  shows  that  C02  emissions  increase  33%  by  2025 
(Dtot  =  1.33)  in  the  BS  scenario,  mainly  due  to  the  increase  of  the 
income,  Dact  =  1 .48,  although  the  energy  intensity  is  reduced  to 
Dint  =  0.88.  In  the  SC-2  scenario  the  emissions  increase  72%  with 
respect  to  2010,  mainly  forced  by  the  strong  increase  of  the  GDP, 
Dact  =  1  -93,  which  is  partially  compensated  by  the  reduction  of 
the  energy  intensity  term,  Dint  =  0.88.  SC-3  scenario  supposes  an 
improvement  with  respect  to  SC-2  scenario,  Dtot  =  1.57,  due  to  the 
reduction  of  the  mixing  term  in  a  11%.  Finally,  SC-4  scenario  has  an 
increase  of  C02  emissions  even  smaller  than  the  BS  scenario, 
Dtot  =  1.31,  which  origin  is  the  small  value  of  the  energy  intensity, 
Dint  =  0.74.  The  coefficients  of  the  LDMI  analysis  are  summarized 
in  Table  4. 

It  is  important  to  note  that  the  reduction  in  the  global  CO 2jnt  is 
twofold,  on  one  hand,  it  is  due  to  the  use  of  a  more  efficient  fossil 
fuel  technology  (lower  energy  intensity)  and,  on  the  other  hand, 
due  to  the  reduction  of  the  share  of  fossil  energy  to  the  energy 
matrix.  Both  contributions  are  equally  important.  The  period 
2011-2025  presents  a  reduction  of  the  global  CO 2,nt  from  2.4  to 
2.1  ktC02/ktoe  in  SC-3  and  SC-4  scenarios. 

5.5.  EKC  verification 

C02  emissions  and  GDP  are  strongly  correlated  in  most  of  cases 
and  the  relationship  between  both  variables  is  expected  to  be 
monotonously  increasing.  The  EKC  suggests  that  the  almost 
linearly  increasing  relationship  between  C02  emissions  and  GDP 
can  be  broken  under  certain  circumstances  as  was  explained  in 
Section  1.  To  face  this  problem  we  first  depict  C02  vs.  GDP,  as  it  is 
shown  in  Fig.  6.  In  this  figure  we  plot  the  raw  and  the  smoothed 
data  for  the  period  1980-2010,  as  well  as  the  results  for  the 
extrapolated  period,  2011-2025,  under  the  four  considered  sce¬ 
narios.  The  black  full  line  that  corresponds  to  the  BS  scenario 
shows  a  moderated  raising  with  a  steadily  increase  of  the  GDP, 
showing  a  clear  linear  tendency  for  GDP  larger  than  250  billion 
USD.  The  red  dashed  curve  that  stands  for  the  SC-2  scenario 
presents  a  more  rapid  growth  of  the  GDP  and  once  more  a  linear 
relationship  with  the  C02  is  fulfilled.  Finally,  green  and  blue  lines 
correspond  to  the  SC-3  and  SC-4  scenarios,  respectively,  and  both 
show  a  rapid  increase  of  the  GDP,  even  faster  than  in  the  SC-2 
scenario,  due  to  the  feedback  mechanism,  but  the  slope  of  the 
curves  is  smaller  than  in  the  BS  and  SC-2  scenarios.  This  suggests 
that  the  latest  scenarios  can  be  on  the  way  to  achieve  a  more 
sustainable  sounding  situation.  In  the  SC-2  scenario  the  GDP 
growths  substantially  and  no  attenuation  measures  are  carried 
out,  that  is  why  the  slope  is  the  same  than  in  the  BS  scenario. 
In  the  SC-3  scenario  the  use  of  renewable  energy  is  promoted  and 
in  the  SC-4  scenario,  moreover,  the  energy  intensity  is  reduced, 
explaining  the  change  in  slope  with  respect  to  BS  and  SC-2 
scenarios.  We  also  mark  in  Fig.  6  the  year  of  the  turning  point 
into  the  second  stage  of  the  EKC  for  the  different  scenarios,  as 
explained  below  in  this  section. 

Prior  to  a  more  quantitative  analysis  it  is  compulsory  to  test  the 
order  of  integration  of  both  variables  LGDPt  and  LC02c  using  the 
tests  of  Ng  and  Perron  [57].7  The  results  in  Table  5  show  that  the 
null  hypothesis  of  no  stationarity  cannot  be  rejected,  indepen¬ 
dently  of  the  statistic  used,  for  both  series,  LGDPt  and  LC02t. 
Accordingly,  both  series  would  be  concluded  to  be  1(1). 


6  Note  that  there  is  an  extra  factor,  Demf,  the  emission  factor,  which  remains 
unchanged  and  equal  to  one. 

7  These  authors  proposed  using  test  statistics,  which  are  modified  versions  of 
Philfip-Perron  and  ADF  tests.  Such  modifications  improve  the  tests  with  regard  to 
both  size  distortions  and  power. 
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Table  4  Table  6 

Results  of  the  C02  emission  decomposition  factors  for  the  period  2010-2025.  Stock-Watson-Shin's  DOLSabcd  estimation  of  linear  cointegration. 


Scenario 

Dfot 

Dact 

Dstr 

Djnt 

Dmix 

BS 

1.33 

1.48 

1.02 

0.88 

0.99 

SC-2 

1.72 

1.93 

1.02 

0.87 

0.99 

SC-3 

1.57 

2.02 

1.02 

0.87 

0.88 

SC-4 

1.31 

1.99 

1.02 

0.74 

0.88 
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Fig.  6.  Venezuela  GDP  vs.  C02  emissions  for  the  period  1980-2025.  Dots  corre¬ 
spond  to  the  official  data,  RF  stands  for  Reference  year  and  TP  for  Turning  point.  (For 
interpretation  of  the  references  to  color  in  this  figure  caption,  the  reader  is  referred 
to  the  web  version  of  this  article.) 
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Table  5 

Ng-Perron  tests'3  for  unit  root. 


MZ™ 

MZ?“ 

MSBols 

mptgls 

Variable 

LGDPt 

-1.458 

-0.497 

0.341 

30.493 

LC02t 

0.361 

0.297 

0.823 

144.345 

Critical  values 

10% 

-14.2 

0.185 

-2.62 

6.67 

5% 

-17.3 

0.168 

-2.91 

5.48 

1% 

-23.8 

0.143 

-3.42 

4.03 

a  The  Modified  Akaike  Information  Criteria  (MAIC)  are  used  to  select  the 
autoregressive  truncation  lag,  k,  as  proposed  in  Perron  and  Ng  [56].  The  critical 
values  are  taken  from  Ng  and  Perron  [57,  Table  1]. 


Once  the  order  of  integration  of  the  series  is  analyzed,  we  will 
analyze  the  long-run  regression  model  [5]  using  the  Dynamic 
Ordinary  Least  Squares  (DOLS)8  estimation  method  of  Stock  and 
Watson  [4]  following  the  methodology  proposed  by  Shin  [58].9 
This  approach  is  similar  to  the  KPSS10 *  tests,  which  are  implemen¬ 
ted  in  two  stages  for  the  case  of  cointegration. 


8  Least  squares  estimation  of  equation  might  suffer  two  problems:  endogeneity 
bias  in  the  explanatory  variables  and  nuisance  parameter  dependencies  due  to 
serial  correlation  in  the  residuals. 

9  In  order  to  overcome  the  problem  of  the  low  power  of  the  classical 
cointegration  tests  in  the  presence  of  persistent  roots  in  the  residuals  of  the 
cointegration  regression,  Shin  [58]  suggests  a  new  test  where  the  null  hypothesis  is 
that  of  cointegration. 

10  These  tests  are  called  Kwiatkowski  et  al.  [59]  tests,  and  assume  the  null 

hypothesis  of  stationarity. 


Parameter 

estimates 

Data  1980- 

2010 

BS  1980- 

2025 

SC-2  1980- 

2025 

SC-3  1980- 

2025 

SC-4  1980- 

2025 

Po 

5.506*** 

.515“* 

7.095*“ 

7.799*** 

8.878*** 

(0.023) 

(0.1777) 

(0.247) 

(0.366) 

(0.513) 

P\ 

1.178*** 

0.764*** 

0.862*“ 

0.723*** 

0.509*** 
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(0.038) 
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R2 
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0.962 

Test:  c; 

0.152 

0.144 

0.187 

0.194 

0.193 

a1 

0.011 

0.028 

0.027 

0.040 

0.057 

*Significance  at  the  10%  level. 

^Significance  at  the  5%  level. 

***  Significance  at  the  1%  level. 

a  Standard  errors  (in  brackets)  are  adjusted  for  long-run  variance.  The  long-run 
variance  of  the  cointegrating  regression  residual  is  estimated  using  the  Barlett 
window  which  is  approximately  equal  to  INT(T1/2)  as  proposed  by  Newey  and  West 

[60]. 

b  We  choose  q  =  INT(T1/3)  as  proposed  by  Stock  and  Watson  [4]. 
c  CA  is  a  LM  statistic  for  cointegration  using  the  DOLS  residuals  from  determi¬ 
nistic  cointegration,  as  proposed  by  Shin  [58]. 

d  The  critical  values  are  taken  from  Shin  [58,  Table  1,  m  =  1]:  =  0.231  (10%), 

C„  =  0.314  (5%),  =  0.533  (1%). 


year 

Fig.  7.  Evolution  of  Venezuelan  C02-GDP  elasticity  for  the  period  2000-2025. 


We  will  estimate  the  coefficients  of  a  long-run  dynamic 
equation  [5]  including  leads  and  lags  of  the  explanatory  variables 
(GDP  in  our  case)  in  the  long-run  regression  model,  i.e.  the  so- 
called  DOLS  regression: 


LC02t=  fiQ+iiALGDPt+  2  fijALGDPt_j  +  6j.  (17) 

j  —  Q 

Moreover  we  use  the  statistic  C^,11  a  LM-type  test  designed  by 
Shin  [58],  to  test  the  null  hypothesis  of  cointegration  against  the 
alternative  of  no  cointegration  in  a  DOLS  regression.  In  Table  6  we 
report  the  estimates  from  the  DOLS  regression  and  the  results 
from  Shin  [58]  test.  Results  show  evidence  of  linear  cointegration 
between  C02  emissions  and  GDP,  because  we  cannot  reject  the 
null  hypothesis  of  cointegration.  The  estimated  value  of  the  C02- 
GDP  elasticity  for  the  dataset  period  is  =1.178,  which  denotes 


11  Cp,  is  the  test  statistic  for  deterministic  cointegration,  i.e.,  when  no  trend  is 
present  in  the  regression. 
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little  environmental  responsibility,  i.e.,  Venezuela  in  the  year  2010 
is  still  in  the  first  stage  of  the  EKC. 

Our  main  goal  is  to  study  the  EKC  of  Venezuela  as  a  function  of 
the  time,  in  the  medium  term,  up  to  2025.  To  carry  out  this  study, 
we  perform  the  same  process  described  above  using  the  projected 
series  provided  by  our  model.  Results  are  shown  in  Table  6.  As  it 
can  be  seen,  for  all  the  scenarios  /i,  <  1  which  proves  that 
Venezuela  is  on  the  way  to  enter  in  the  second  stage  of  the  EKC. 
In  the  case  of  the  BS  scenario,  the  elasticity  goes  from  =  1.178  in 
2010  to  0.764  in  2025,  while  in  the  SC-2,  SC-3,  and  SC-4  scenarios  it 
moves  till  0.862,  0.723,  and  0.509  in  2025,  respectively.  In  Fig.  7  we 
show  the  evolution  of  the  C02-GDP  elasticity  as  a  function  of  the 
year  for  the  period  2000-2025.  This  figure  shows  a  very  stable 
value  of  the  elasticity  for  the  reference  period,  up  to  2010,  and 
then  a  down  sloping  tendency  in  all  the  considered  scenarios, 
which  means  that  Venezuela  is  moving  into  the  second  stage  of 
the  EKC.  In  the  BS  scenario  the  elasticity  goes  rapidly  below  1 
(turning  point  in  2014)  but  remains  stable  for  the  rest  of  the 
period.  To  understand  this  fact  we  have  to  take  into  account  that 
the  data  concerning  energy  intensity  in  Venezuela,  which  shows  a 
quite  erratic  evolution  in  the  reference  period  in  all  the  industrial 
sectors,  leads  to  a  down  sloping  smooth  projection  for  the  period 
2011-2025.  This  reduction  of  the  energy  intensity  is  the  origin  of 
the  reduction  of  the  elasticity.  The  same  argument  holds  in  the  SC- 
2  scenario,  because  for  this  scenario  the  energy  intensity  has  the 
same  evolution  than  in  the  BS  scenario.  In  this  case  the  entrance  in 
the  second  phase  of  the  EKC  happens  latter,  turning  point  in  2016. 
The  behavior  of  the  SC-3  scenario  is  very  similar  to  SC-2  scenario, 
although  in  this  case  renewable  energy  is  promoted.  The  turning 
point  of  the  SC-3  scenario  also  happens  in  2016.  Finally  in  the  SC-4 
the  faster  decreasing  of  the  elasticity  is  obtained,  showing  the 
turning  point  in  2014.  Besides  having  the  different  scenarios 
different  turning  points,  they  also  have  a  different  asymptotic 
behavior.  BS  and  SC-2  scenarios  seem  to  converge  to  a  constant 
value,  around  /./,  =  0.76,  while  in  SC-3  and  SC-4  scenarios  the 
elasticity  continues  decreasing  all  the  time.  We  can  extract  as  a 
conclusion  that  although  it  is  possible  to  reduce  the  rate  of 
emissions  even  without  implementing  any  sustainable  goal, 
rapidly  the  country  gets  a  saturation  value.  However,  undergoing 
the  appropriate  sustainable  measures  the  country  can  steadily 
reduce  the  C02  emission  rate  and  eventually  reduce  C02  emis¬ 
sions.  The  obtained  results  are  not  general  and  depend  very  much 
on  the  country,  i.e.,  of  the  reference  data  set.  Indeed,  in  a  previous 
work  for  Ecuador  [2]  the  evolution  of  the  value  of  the  elasticity  is 
quite  different  to  Venezuela. 

Changes  made  in  the  SC-4  scenario  (promotion  of  renewable 
energy  and  increase  of  the  energy  efficiency  in  the  productive 
sectors)  help  to  create  a  more  favorable  future  scenario  in  terms  of 
sustainable  development.  Undergoing  the  goals  of  this  scenario, 
the  country  could  be  able  to  enter  the  area  of  environmental 
stability  in  the  medium  term. 


6.  Summary,  conclusions  and  policy  implications 

In  this  paper  we  have  estimated  the  value  of  the  GDP,  the 
energy  consumption  and  the  C02  emissions  in  Venezuela  in  the 
medium  term,  from  2011  to  2025.  The  model  was  based  on  a  nexus 
relationship  which  is  a  variation  of  the  Kaya  identity  and  on  a  GDP 
formation  which  includes  the  effect  of  renewable  energy.  This 
model  was  first  presented  by  Robalino-Lopez  et  al.  [1,2]  and 
applied  to  Ecuador.  The  official  data  set  for  the  period  1980- 
2010  was  used  to  calibrate  the  model,  while  in  the  period  2011- 
2025,  an  estimation  of  different  variables,  including  C02  emissions 
and  sectoral  energy  consumption,  was  carried  out.  On  the  other 
hand,  different  scenarios  that  represent  the  evolution  of  GDP, 


energy  matrix  and  productive  sectoral  structure  have  been 
defined.  One  of  the  key  points  of  the  model  is  the  direct  influence 
of  renewable  energies  on  the  GDP  [18],  which  creates  a  feedback 
mechanism  between  income  and  C02  emissions.  The  relationship 
between  C02  emissions  and  GDP  is  studied  on  the  light  of  the 
environmental  Kuznets  curve  (EKC)  and  the  value  of  the  C02-GDP 
elasticity  is  calculated  as  a  function  of  the  time.  It  is  found  that 
Venezuela  will  enter  into  the  second  phase  of  the  EKC  in  the 
coming  years. 

We  have  carried  out  the  projection  of  C02  emissions  under  four 
different  scenarios.  First,  in  the  BS  scenario  (baseline  scenario),  the 
variables  of  the  model  were  taken  accordingly  to  the  observed 
tendency  in  the  reference  period  1980-2010.  For  this  scenario  the 
emissions  in  2025  are  33%  higher  than  in  2010.  Similarly,  the 
second  scenario,  SC-2,  is  characterized  by  the  doubling  of  income 
in  2025  relative  to  2010.  In  this  case  the  emissions  are  72%  higher 
than  in  2010.  In  the  third  scenario,  SC-3,  besides  assuming  the 
doubling  of  the  GDP,  we  impose  a  decrease  of  the  fossil  energy 
share  down  to  76%,  reaching,  C02  emissions  at  the  end  of  the 
period  57%  higher  than  in  2010.  Finally,  in  the  fourth  scenario,  SC- 
4,  we  also  include  changes  in  the  productive  sectoral  structure  to 
achieve  a  reduction  of  energy  intensity,  which  supposes  a  decreas¬ 
ing  of  the  C02  intensity.  In  this  case,  the  emissions  are  even  lower 
than  in  the  BS  scenario,  only  31%  higher  than  in  2010.  In  the  SC-3 
scenario  the  increasing  use  of  renewable  energy  allows  to  reduce 
C02  emissions  15%  with  respect  to  SC-2  scenario.  Finally,  in  the  SC- 
4  scenario  an  improvement  of  the  energy  intensity  is  implemen¬ 
ted,  reaching  a  reduction  of  41%  with  respect  to  the  SC-2  scenario. 
Therefore,  we  can  conclude  that  the  reduction  of  C02  emissions 
can  be  accomplished  using  two  different  strategies:  improving  the 
energy  intensity  and  increasing  the  share  of  renewable  energies  to 
the  total  energy  consumption.  Because  of  the  lack  in  the  literature 
of  projections  of  C02  emissions  for  Venezuela,  we  cannot  directly 
compare  our  results  with  any  other,  although  one  can  find  in  the 
bibliography  some  reference  points  with  which  to  check  the  model 
and  the  scenarios.  On  one  hand,  our  estimations  of  C02  emissions 
are  larger  than  the  target  proposed  by  the  National  Electricity 
Corporation  of  Venezuela  for  2050,  which  corresponds  to  half  the 
emissions  in  1999  [41],  This  target  is  based  on  the  renewable 
energy  resources  available  in  Venezuela  and  not  yet  used.  There¬ 
fore  the  C02  emissions  for  the  SC-3  and  SC-4  scenarios  could  be 
even  lower  if  the  share  of  renewable  energy  to  the  energy  matrix 
is  increased  more  than  in  the  proposed  scenarios.  On  the  other 
hand,  in  [23]  it  is  shown  that  the  C02  efficiency  of  Venezuela 
increased  from  1997  to  2007,  which  suggests  that  SC-3  and  SC-4 
are  probably  rather  realistic  scenarios. 

Concerning  the  analysis  of  the  EKC  hypothesis  we  conclude  that 
Venezuela  do  not  fulfill  this  hypothesis  under  any  of  the  consid¬ 
ered  scenarios,  i.e.,  in  the  medium  term  the  C02-GDP  elasticity  is 
always  positive.  Our  results  show  that  Venezuela  in  2010  is  still  in 
the  first  stage  of  the  EKC  and  that  under  the  BS  scenario  can  only 
slightly  reduce  the  value  of  the  elasticity  and  it  will  remain  with  a 
rather  stable  value  if  no  attenuation  measures  are  implemented  in 
the  near  future.  Very  much  the  same  situation  happens  under  the 
SC-2  scenario,  which  only  differs  of  the  BS  scenario  in  the  rapid 
increase  of  the  GDP.  Although  the  country  enters  into  second  stage 
of  the  EKC  relatively  soon,  in  2014  and  2016  for  the  BS  and  SC-2 
scenarios,  respectively,  the  elasticity  moves  into  a  constant  value, 
/tj  rs  0.76,  in  the  medium  term  and  therefore  the  country  will 
move  towards  a  dead  end,  concerning  environmental  sustainabil¬ 
ity.  In  SC-3  and  SC-4  scenarios  different  environmentally  sounding 
measures  are  implemented  and,  once  more,  the  country  enters 
very  rapidly  into  second  stage  of  the  EKC,  also  in  the  years  2016 
and  2014  for  the  SC-3  and  SC-4  scenarios,  respectively.  However,  in 
this  case  the  trend  of  the  elasticity  is  clearly  down  sloping  and  it 
could  be  expected  to  reach  a  negative  value  in  the  medium  term, 
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especially  in  the  SC-4  scenario.  The  diversification  of  energy 
sources  and  the  reduction  of  the  energy  intensity  are  essential  to 
reach  a  sustainable  development.  Fossil  fuel  substitution  has  been 
proven  to  be  a  successful  way  to  reduce  emissions  by  moving 
energy  consumption  into  forms  with  lower  carbon  intensity  and 
improving  energy  efficiency.  Our  results  are  in  agreement  with 
Narayan  and  Narayan  [25]  who  found  that  Venezuela  for  the 
period  1980-2004  showed  evidences  of  the  fulfillment  of  the  EKC, 
with  a  long-run  elasticity  of  0.29,  much  smaller  that  the  short-run 
one,  1.57,  which  is  a  hint  of  the  accomplishment  of  the  EKC 
hypothesis. 

The  methodology  presented  in  this  paper  can  be  useful  to 
estimate  C02  emissions  of  a  given  country  and  to  understand  the 
driving  forces  that  guide  this  process,  such  as  economic  growth, 
energy  use,  energy  mix  structure,  and  fuel  used  in  the  productive 
sectors.  Moreover,  it  can  serve  as  a  pedagogical  tool  for  explaining 
to  policymakers  the  possible  ways  to  design  a  policy  for  reducing 
C02  emissions  in  a  medium  term  horizon.  As  was  stated  by 
Robalino-Lopez  et  al.  [1,2],  the  presented  methodology  is  transfer¬ 
able  to  other  period  of  time  or  region. 
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